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t— h Abstract 

An incomplete sky coverage poses difficulties in the extraction of the weak lensing information 
from the CMB. We test the reconstruction of the weak lensing convergence from CMB maps to 



o 

which masks of point sources have been applied. We use the quadratic estimator with a kernel 
with finite support acting in real space for a Planck simulation. We recover the lensing signal 

p 

without significant loss of power or addition of spurious correlations, thus showing that masking 
defected pixels does not affect the reconstruction of the weak lensing convergence in real space. 
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I. INTRODUCTION 



The weak gravitational lensing of the CMB is a measure of the gravitational potential 
projected along the line of sight from the last scattering surface. It was indirectly detected 
using WMAP data [H [2] and recently direct detections in WMAP [3| and ACT [4] were 
reported. 

In a CMB experiment, the temperature signal is mixed with various foreground signals 
which can generate spurious correlations. To obtain the signal of interest, we must first 
remove the contaminating signals. However, an incorrect removal could create new correla- 
tions. In order to robustly test estimators of quantities measured from the CMB, we need 
to ensure a correct removal of the contaminants without degrading the signal of interest. 

Any process resulting from the coupling of different Fourier modes is sensitive to contam- 
inants. The weak gravitational lensing of the CMB induces non-Gaussianities in the form of 
mode couplings and consequently is very sensitive to spurious correlations from foreground 
signals. The contaminants can be either physical components (integrate Sachs- Wolfe effect, 
Sunyaev-Zel'dovich effect, galactic foregrounds, point sources) or experimental components 
(instrumental noise, beam effects). Such systematics have become more relevant with the 
current generation of high-precision CMB experiments such as ACT 1 , Planck 2 or SPT 3 . 
Here we are concerned with the treatment of point sources on CMB maps. 

Point sources can be removed either by subtracting the estimated emission, through a 
model of the point source flux, or by discarding the contaminated pixels above a certain flux. 
Filtering the emission from point sources was the method used in the recent detection of the 
convergence power spectrum from the ACT data jl]. Another approach is called inpainting 
and consists in removing the point sources and replacing them with an extrapolation of the 
missing data using some prior of the complete data j5j [6] . However, both the subtraction 
of the point source emission and the inpainting of masking disks are dependent on the 
extrapolation model (i.e. the model of the point source flux and the prior of the complete 
sky, respectively). To extract a model independent result, we should be able to simply 
discard the contaminated pixels by applying a point source mask to the map we want to 
analyse. 

Discarding the pixels amounts to multiplying the map by a mask, which causes spatial 
modulation of the average sky brightness and consequently affects the analysis of CMB 
experiments. The multiplication of the map by a mask for all point sources can be re- 
garded as the successive multiplication of the map by a mask for each point source [?]. 
The Fourier transform of the masked map will be equivalent to successively convolving the 
Fourier transform of the original map by that of a series of such masks, which will mix power 



1 |http : / /www . physics .princeton. edu/act/about .html 

2 http : //www.rssd. esa. int/index .php?project=planck 

3 http : //pole . uchicago . edu/ 
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of the map over a range of multipoles. The reconstructed power spectrum must then be 
carefully corrected from these contaminants introduced by the analysis of incomplete data. 

Here we want to study to which extent masking small areas at the locations of point 
sources affects the weak lensing convergence extracted from CMB maps with the estimator 
introduced in Ref. [8]. This is a variant of the quadratic estimator with the novelty of having 
a finite support and acting in real space. (See also Ref. [9j for a study of its properties.) 
We also compare the results with those obtained from the conventional quadratic estimator 
defined in harmonic space. 

The manuscript is organized as follows. In section |TT] we describe how we estimate the 
lensing convergence using the minimum variance estimator and how we deal with the bias of 
the estimator. In section III we synthesise the point source masks and apply the estimator 
to lensed CMB maps with and without masking. In section IV we study the degradation 
of the estimation of the convergence power spectrum as a function of the parameters of the 
mask, namely the point source number density and the masking disk diameter per point 
source. We summarize our results in section [V] 



II. ESTIMATOR OF THE WEAK LENSING CONVERGENCE 

Weak lensing of the CMB consists in the deflection of CMB photons from the original 
propagation direction 6 on the last scattering surface (the source plane) to an observed 
direction on the sky today (the image plane) (see Ref. pU] for a review). This deflec- 
tion amounts to remapping the unlensed temperature anisotropics T to the lensed ones T 
according to 

f{6) = T{0) = T(0 + W) = T{6) + W> ■ VT(0) + 0[(W) 2 ]. (1) 

The deflection angle ex = — is given by a = Vip, where the lensing potential ip is the 
projection, on the image plane and along the line-of-sight, of the gravitational potential of 
the intervening large-scale structure, and V is the covariant derivative on the image plane. 
The right-hand-side of Eq. Q is written for small deflections as a first-order Taylor expansion 
in the deflection, showing that the lensed CMB temperature depends on the derivative of 
the unlensed temperature through the coupling with the deflection. 

For small deflections, the lensing effect is described by the amplification matrix A = 
1 — k, where the convergence tensor n is defined by second-order derivatives of the projected 
potential. In the absence of rotations, k, can be decomposed as a sum of an isotropic 
(diagonal) and an anisotropic (traceless and symmetric) part as follows 

/ K + K + K x 

K = 

V K x Kq — K + 

The isotropic part produces a convergence/ dilation whereas the anisotropic part produces 
a shear. In particular, the convergence k = — (d 2 . + dy)ip/2 magnifies a feature on the 
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last-scattering surface, the shear k + = — {d 2 — dy)ip/2 stretches it along the x-axis while 
compressing it along the y-axis, and analogously the shear k x = —d x d y ip stretches it along 
the y = x axis and compresses it along the y = — x axis. 

From Eq. (fl|) the lensed temperature power spectrum writes 

(f(£') f(£ -£')) = (2tt) 2 6 b {£) Op + (2tt) 2 [£ ■£' C e , + £■(£- £') C ]t _ t >\] i/>(£). (3) 

The lensing field can thus be reconstructed with a quadratic estimator [TT], i.e. a weighted 
convolution of the square of the lensed map T{£). The lensing information is contained in 
the off-diagonal terms of the lensed T(£)T(£') correlation that are generated by anisotropies 
induced by the lensing potential. In Ref. [T2] an optimal weighting is derived as well as a 
normalization that, in the assumption of Gaussian temperature and lensing random fields, 
and Gaussian noise uncorrelated with the signal, ensures the estimator has no calibration 
bias. 

Here we will use the real-space estimator introduced in Ref. [8] , which is a real-space ana- 
log of the optimal estimator based on the convolution of the square of the lensed temperature 
map with a local kernel as follows 

k (o) = j d 2 e' f(o') J d 2 e" f(o") Q{e,e',e"). (4) 

Here the kernel Q(6,6' ,6") is related to that of the standard quadratic estimator Q(£,£') 
by 

Q(6,e',e") = J ^0 vq>[i£-0] J 0^ e^[-i£' ■e']e^[-i{£-£') ■ 0"} Q(£,£'). (5) 

The kernel Q(£, £') is the minimum variance (to leading order) weight function that convolves 
the temperature map in the estimation of the weak lensing convergence in harmonic space 

ko(i)=i 2 j 0r 2 f{£') f{£-£!) Q(£,£'). (6) 

For comparison of the results, we will also use the harmonic-space estimator. The estimators 
will be applied to lensed CMB maps synthesised in the flat-sky approximation. 



A. The estimated map 

The reconstructed convergence maps contain not only the contribution due to the lensing 
potential but also a contribution due to the unlensed CMB temperature, and we write it 
schematically as 

— K o|v> + K o\tp=o- (7) 
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The non-lensing contribution is due to other anisotropic couplings in the CMB tempera- 
ture field, such as T VT correlations. These vanish when taking the average over realizations 
of the CMB (i.e. realizations of lensed CMB maps with the same lensing potential, being 
thus a source of noise), but are not zero in a given realization. A convergence map recon- 
structed from one unlensed CMB realization has a rms per pixel of cr reo|i0=o = 0.1. This is 
the same amplitude as the rms of a convergence map reconstructed from a lensed CMB 
map, showing that the lensing contribution is sub-dominant. Over 100 realizations, the rms 
reduces to & Kol4l=0 = 0.01. Hence the estimator is noisy but unbiased. 

However, the estimator is biased if some assumptions are dropped. Firstly, if all config- 
urations of £, £' and {£ — £') are used when computing Eq. the diagonal terms of the 
lensed T(£)T(f) correlation will also contribute to Ko|v>=o by introducing an additive bias in 
the estimator. Since both jcqu an d the noise should have zero mean when averaged over the 
map, we may subtract the pixel average, k — (kq) - el , to remove that zero-point function. 
Secondly, the estimator reconstructs the lensing contribution kq^ with a calibration bias 
if non-Gaussian structures contribute to the lensing field |13j . Although the CMB lensing 
efficiency peaks at redshift z ~ 3, where most structure is linear, it has contributions from 
all redshifts so that inevitably it will be sourced by non-Gaussian structures. 

The reconstructed maps depend on the CMB detector noise via the kernel in the estimator 
(see Ref. JHj). The power spectrum of the CMB temperature detector noise depends on the 
characteristics of the simulated CMB experiment. For one frequency channel, the CMB 
noise power spectrum is given by [HI US] 

N(£) = (6 {whm a pix ) 2 expi(e iwhm ) 2 £(e+ l)/(81n2)]. (8) 

Here 6{ w hm is the beam full-width at half-maximum, which defines the pixel size in the 
simulated maps, and <7 P i x is the white noise amplitude per pixel in units /zK, which depends 
on the detector sensitivity and observation time. Our maps simulate a fraction of the sky 
/ s ky smaller than the fraction of the sky covered by the target experiment f^ y . In order 
to have the same noise amplitude in the maps, we adopt the so-called active incomplete 
sky-coverage strategy [16], which assumes the smaller coverage is compensated by longer 

1/2 

observation times, reducing the noise per pixel as <T pix = (/sk y // s ky) • Alternatively, 
this can be seen as scaling the pixel size such as to cover the fiducial fraction of the sky. 

B. The estimated power spectrum 

We turn now to the power spectrum of the convergence map. The noise in the map is 
assumed to be uncorrelated with the lensing signal, (/«o|^ K o|v>=o) = ^> but it has a non- 
vanishing power spectrum (k \^ =0 «qi^_ ). Consequently, the convergence power spectrum 
computed from an unbiased convergence map is biased. This additive bias is evaluated 
in Ref. [IT] and can be subtracted off. Alternatively, we estimate convergence maps from 
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Figure 1: Error budget of the lensing power spectrum estimation for a PLANCK-like 



survey. The three terms in the square-root of Eq. (11) are shown in dashed-dotted (estimator 
noise), solid (cosmic variance), and dashed (bias uncertainty) lines. 



various realizations of unlensed CMB maps for the same cosmology. The average of the 
power spectra over the different CMB realizations is an estimator of (ko|^=o ^q| v=o) • ^ a ^ so 
includes the additional bias in the case of a non-zero zero-point. We thus subtract it off of 
the estimated biased power spectrum to obtain a corrected estimator of the power spectrum 
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However, this estimator is still biased. Indeed, since C^° K ° is a four-point correlator of 
the lensed CMB field, it also includes a connected part which introduces couplings between 
lensing modes, thus producing an additive bias which depends on the convergence 
power spectrum itself [17]. Further additive biases arise when expanding Eq. ([I]) to higher 
orders in the lensing potential [18J. The magnitude of these higher-order biases was studied 
numerically in Ref. [T3] . 

Here we will use Eq. (|9| to compute the convergence power spectrum, thus correcting for 
the lensing-independent noise-bias, and will also check for a multiplicative bias. We include 
three contributions to the variance of the corrected power spectrum, arising from: the 
variance of the convergence estimator (the noise), the cosmic variance and the uncertainty 
on the removed bias. The theoretical variance of the estimator per binned mode Var[K ] is 
derived in Ref. [8], Appendix A. Adding the cosmic variance contribution yields the variance 
of the biased estimator of the power spectrum. The uncertainty on the additive bias is the 
variance of the mean (C^ 0| ' / ' =oKo, ' / ' =0 ) CMB , which was estimated by averaging the convergence 
power spectra over N realizations of unlensed maps. It is given by 



Var[C, 
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where, for each map, the variance Var[C^ 0| ' / ' =oKo| ' / ' =0 ]j is computed as in the first term of 
Eq. ( 11 ). The variance of the corrected power spectrum is obtained by adding the contribu- 
tions in quadrature, written in terms of variance per band as follows 

VarfCf 01 *' 801 *] « j-j— (Var[£ ] + Cf^) 2 + Var[C; o| *=° K(W=0 ]. (11) 

" /sky 

The denominator accounts for the number of modes probing the multipole band, which 
depends on the area of the simulated maps. To forecast the error budget of a given experi- 
ment, its fiducial fraction of the sky must be used instead. The variance-covariance matrix 
is essentially diagonal [T7j, i.e. the multipoles may be considered independent. 

The three contributions are shown in Fig. [T] for a PLANCK-like survey with the noise 
specifications of Tab. [I] and / s k y = 1. The noise term dominates on all scales. The bias is 
removed using 100 realizations, which is a large-enough sample to produce a sub-dominant 
bias uncertainty. 



III. APPLICATION TO MASKED LENSED MAPS 
A. The point source mask 

We generate the map 

M = [MB + N] [SB] , (12) 

where M, B, N and S refer respectively to lensed CMB, beam, noise and point source maps, 
as described below. 

We start by generating a lensed CMB map using the formalism described in Ref. [8], 
Appendix B. In summary, we first compute the TT and the ipip power spectra generated 
by CAMB [19], on scales I < 4000, for a fiducial flat ACDM cosmology consisting of h = 
0.7, tt cdm h 2 = 0.112, tt h h 2 = 0.0226, n s = 0.96, a 8 = 0.85, r rci = 0.09, T = 2.725 K. With 
these power spectra, we build a CMB temperature map and a deflection map on a square 
patch. We then apply the deflection map to the CMB map to produce a lensed CMB 
map. This operation consists in shifting the CMB map by the lensing potential assuming 
the Born approximation. We then introduce the detector effects, first by convolving the 
lensed map with the beam profile and then by adding a map constructed from the detector 
noise power spectrum. The size of a map is given by the pixel size times the number of 
pixels, with the pixel size being determined by the beam size. We use two cases based on 
the experimental characteristics of Planck, denoted by Designer and Planck, the first 
containing no detector noise, the second containing that of the v = 143 GHz Planck 
channel of the HFI instrument [20]. Finally, we mask the map by placing disks of finite 
diameter at the point source locations, where the amplitude is set to zero. 
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pix 


rM 

sky 


FOV map 


n s 


0s 




(GHz) 


(arcmin) 


(uK) 




(deg 2 ) 


(deg" 2 ) 


(arcmin) 


Designer 


143 


7.2 





1 


61.3x61.3 


0.03 


36 


Planck 


143 


7.2 


7.5 


1 


61.3x61.3 


0.03 


36 



Table I: The specifications of the two experiments used in this study. The experiment 
denoted by Planck is inspired by the v = 143 GHz channel of the Planck HFI instrument |20j . 
The experiment denoted by Designer is the same as Planck for no detector noise. 

The mask is generated by placing a number density n s of excisions on locations of the 
map randomly generated with a uniform distribution. The value of n s used for the cases of 
Planck and Designer is obtained from the total number of point sources detected on the full 
sky at the 143 GHz channel, as reported in the Planck early results [2T]. The mask should 
be convolved with the beam, which is equivalent to assigning to each point source a disk of 
finite diameter 9 S . The size of the masking disk 9 S , assumed uniform, must be proportional to 
the size of the beam ^fwhm- Indeed, a higher-resolution experiment, i.e. one with a smaller 
beam, would smear the flux of each source over a smaller area and in addition it would 
detect more point sources. Hence, if the sources are uniformly distributed, the product 
n s 6 s is constant across experiments. For example, the analysis of the Atacama Cosmology 
Project (ACT) 2008 Survey used masks with 8 S = 5 6>f wmT i [22] . for the n s = 0.35 deg -2 point 
sources detected with a beam 6*fwhm ~ 2 arcmin [23] . This implies masking disks of the same 
relative size, 8 S = 5 #f w h m > for Planck . Thus, in our masks we will use for the masking disk 
diameter the value 8 S = 36 arcmin. We note that pixels where masks overlap count as one 
point source only. 

B. The convergence of a masked lensed map 

We apply the real-space estimator to a lensed CMB map with and without the point 
source masks, and obtain estimated maps of the convergence which we denote respectively 
by Rq and K o (whenever there is no case for confusion, to simplify the notation we drop 
the hat that identifies an estimated quantity). For comparison, we also implement the 
harmonic-space estimator on the same maps (with the implementation also described in 
[8]). In addition, we run several realizations of unlensed CMB maps. Each realization 
consists of a different map built from the same temperature power spectrum. For each 
unlensed CMB map, we apply various realizations of the mask with the same n s and 8 S but 
with the excisions located at different positions. 

The reconstruction made with the real-space estimator for a masked map, ko, is very 
similar to the reconstruction kq, only missing a small amount of the structure obtained in 
k . The difference between the two maps k — k (shown in Fig. [2j left panel) is confined to 
localized medium-scale structure similar in scale and location to that of the mask used in this 
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Figure 2: Maps for the Planck experiment. Left panel: the difference (kq — kq) between the 
two maps reconstructed with the real space estimator is shown in a greyscale ranging from black 
(no difference) to white (largest difference). Central panel: the mask map, with amplitude on the 
white dots representing the effective size of the point sources, and amplitude 1 elsewhere. Right 
panel: the (kq — kq) map obtained from the harmonic space estimator, same scale as in the left 
panel. 

realization (shown for comparison in Fig. [2j central panel). For the case of the harmonic- 
space estimator the result is very different, as shown in Fig. [2] (right panel). Indeed, due to 
the highly nonlocal character of the harmonic filter, there is a spread of artifacts from the 
disks over the whole map, producing a large amount of new structure which overwhelms the 
lensing structure. This is alternatively seen in the power spectrum of the difference maps, 
shown in Fig. [3] (left panel). The power spectrum of the difference map produced with the 
harmonic-space estimator without detector noise increases at all scales. Adding detector 
noise, this increase of power is predominant at small scales. In contrast, the power spectrum 
of the difference map in the real-space case is peaked, showing an increased difference on 
scales I ~ 200, which corresponds to roughly three times the size of the masking disks, and 
implying a localized small spread of the artifacts. The histogram in Fig. [3] (right panel) 
shows that 97% of the pixels have an identical value in kq and k,q. The remaining 3%, which 
corresponds to three times the area masked according to Tab. [T| are distinctly grouped at 
values lower than the typical amplitude of a noisy pixel a K0 =0.1. 

We note that the comparison of the power spectra shapes in Fig. [3] (left panel) is made 
with arbitrary amplitudes. This is because the normalization in the estimators fails in the 
presence of noise and masking, and the maps become affected by different calibrations. The 
multiplicative bias of a map can be evaluated by taking the ratio of the cross power spectrum 
(jKK in ^ e au ^- p Qwer spectrum of the input convergence map C^ inKin . Assuming that the 
noise in the reconstructed power spectrum is uncorrelated with the lensing information, this 
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Figure 3: The (kq — kq) map. Left panel : The power spectrum of the difference map for the real 
space (black) and harmonic space (gray: dashed-dotted for Designer, full for Planck) estimators, 
with arbitrary amplitude. Right panel : The histogram of the difference map obtained with the 
real space estimator. 



ratio is a direct measure of the calibration bias. 



C. The convergence power spectrum of a masked lensed map 

We compute now the convergence power spectrum and its variance from the convergence 
maps estimated from lensed CMB maps, both in the case of masked and unmasked maps 
and for the two experiments of Tab. [TJ The power spectra are computed for all £-modes in the 
range 6 < i < 1500. The small-scale limit corresponds to twice the pixel size, according to 
the Nyquist sampling theorem, and the large scales are limited by the size of the maps. We 
follow the procedure described in Sec. [TTJ subtracting the mean convergence power spectrum 
from several realizations of the unlensed sky, as described by Eq. ([9]), which succeeds in 
removing most of the additive bias. We then apply a calibration derived from the ratio 
of the estimated power spectrum (after subtracting the additive bias) to the input power 
spectrum. In particular, we use the case Designer as the control experiment for calibration, 
and average this ratio over the region where the two power spectra are approximately parallel 
i.e. 70 < i < 400. After this correction, we remain with the residual multiplicative bias 
shown in Fig. [4j The residual, shown with no binning, fluctuates around 1 in the range 
70 < i < 400, with the fluctuations being larger in the case of the harmonic-space estimator. 

The results for the debiased and calibrated convergence power spectrum are shown in 
Fig.[5]for Designer and in Fig.[6]for Planck. The power spectrum is binned in logarithmically 
spaced intervals and the error bars are correspondingly reduced according to the number of 
£-modes per bin. The error bars include the effective detector noise per pixel, as discussed 
earlier on, and are produced for the effective size of the maps. Forecasted error bars for the 
Planck survey are smaller by a factor of 3 and the power spectrum of the estimator noise 
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Figure 4: The multiplicative bias. The residual multiplicative bias in the power spectra for 
the real-space (black line) and harmonic-space (gray lines) estimators. 
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Figure 5: Estimated convergence power spectrum for Designer. Left panel : Pristine 
lensed map. Right panel : Masked lensed map. In both cases are shown the reconstruction with 

the harmonic-space estimator (gray diamonds), with the real-space estimator (black squares), the 
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for that case (assuming /J^ = 1) is also shown. 

The real-space estimator produces similar results for all the four cases (whether masked 
or unmasked, with or without detector noise). The four power spectra obtained show similar 
shapes and amplitudes and the same calibration is applied to all. In particular, the recon- 
struction on the masked maps is remakably similar to the reconstrucion on the unmasked 
maps, just slightly lower around £ ~ 200 due to loss of power in the masked regions seen 
in Fig. [3] (left panel). Concerning the reconstruction on maps in the presence of detector 
noise, the estimator is insensitive to uncorrelated noise, since the sum of the product of pairs 
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Figure 6: Estimated convergence power spectrum for Planck. Left panel : Pristine lensed 
map. Right panel : Masked lensed map. In both cases are shown the reconstruction with the 
harmonic-space estimator (gray diamonds) , with the real-space estimator (black squares) , the input 
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of neighbouring pixels, weighted by the kernel, averages out uncorrelated pixel noise. Also 
notice that for all cases the recovered power on the smallest scales is limited by the size of 
the kernel (£ = 1000), visible as a decrease in the last two bins. Indeed, as discussed in 
Ref. [8], there is an averaging of the modes smaller than the extent of the kernel. 

The behaviour of the harmonic-space estimator is very different. It performs well only 
in the case of no detector noise and no masking. In the presence of detector noise, it picks 
up the white noise, as it was shown in Ref. [8], increasing the power on the noise-dominated 
small scales. In the presence of masking, the spurious power changes both the shape and 
amplitude of the estimated power spectrum, and the calibration found with the pristine map 
no longer applies. 



IV. VARYING THE MASK PARAMETERS 

We have two attributes to characterize the point source mask: the number density of 
sources and the masking area per source. To understand how these two parameters affect 
the lensing reconstruction, we produced maps where we either a) varied the number density 
of excisions while keeping the diameter of the masking disks constant at 9 S = 5 9{ w hm, or b) 
varied the diameter while keeping the number density constant at n s . 

We start by assessing the effect of the number density of masked sources in 
the estimation of the convergence. We produce maps with S = 5 iwhm for n s G 
{0, 0.01, 0.03, 0.09, 0.18, 0.3, 0.6, 0.9, 1.2, 1.5, 1.8, 2.4, 3.0, 3.3, 4.2, 6.0, 9.0, 12, 15, 21} deg~ 2 , 
which corresponds to iV s G {0, 38, 113, 338, 676, 1126, 2252, 3379, 4505, 5631, 6757, 9009, 11262, 
12388,15767,22524,33785,45047,56309,78833} sources. For each value of the number 
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Figure 7: Degradation of the reconstructed convergence power spectrum with the mask 
parameters. Left panel : The normalized difference between the power spectrum modes I = 103 
reconstructed with and without masking with the real-space estimator for the case of Planck, are 
shown as a function of n s for each of the 10 realizations (gray squares) and averaged over the 10 
realizations (black squares). Right panel : As the left panel, but as a function of 9 S . 
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density, we generate ten realizations of the mask which we then apply to one realization of 
the lensed CMB map. We implement the real-space estimator to the resulting masked maps, 
and present the results in Fig. [7] (left panel) showing the fractional difference between the 
power spectrum estimated with each of the masks and the unmasked case. The difference 
is normalised with the power spectrum estimated for the unmasked case, and is shown for 
the mode £ = 103. The fiducial case (for Planck) is the point at n s = 0.03 deg -2 . 

We observe that the estimation of the convergence power spectrum from a map which 
has been masked degrades monotonically with the increase of the density of excisions in the 
mask as compared to the estimation from the complete map. Both the Planck and Designer 
(not shown) cases show similar behaviours, with the reconstructions in the presence of 
noise detector showing a larger dispersion. The amplitude of the power spectrum in the 
presence of masking decreases with the masking density due to the loss of power discussed 
In particular, in a (x, ln[y]) plot, ln[C^ 103 ] is a negative linear function of the 
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excision number density n s , which implies that the degradation may be written in the form 
C^ 103 (n s ) = C^ 103 (0) exp[— an B ]. The monotonically decrease of ln[C^ 103 ] implies that the 
degradation tends to 100%, which is reached when the mask covers the whole map. 

We proceed to assess the effect of the size of the masking disks in the esti- 
mation of the convergence. We generate point source masks with different mask- 
ing disk diameters 8 S while keeping the number density constant. In particu- 
lar, we use n s = 0.03 deg -2 , which corresponds to 113 sources, for 8 S G 
{0, 1, 2.5, 5, 10, 15, 20, 25, 30, 35, 40, 50, 60, 70, 80, 90, 100, 110, 120, 140} 9 lwhm . For each value 
of the diameter, we generate ten realizations of the point source mask which we then apply 
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Figure 8: Amplitude degradation. Normalized difference between the power spectrum modes 
t = 103 reconstructed with and without masking with the real-space estimator as a function of the 
masked fraction, for varying n s (solid) or varying S (dashed-dotted) . 



to the same realization of the lensed CMB map. We implement the real-space estimator to 
the resulting masked and unmasked maps, and show the resulting fractional deviations in 
Fig. [7] (right panel), where the fiducial case (for Planck) is the point at 9 S = 5. 

As before, we observe that the estimation of the convergence from a map which has been 
masked degrades with the size of the masking disks as compared to the estimation from the 
complete map. In a (x, ln[y]) plot, ln[C^ 103 ] is now a negative quadratic function of the disk 
diameter 9 S , which implies that it can be written as C^ W3 (9 S ) = C^ 103 (0) exp[— 9 2 /(2a 2 )]. 

We can directly compare the two analyses if we translate n s and 9 S into masked fractions. 
Fig. [9] shows this comparison in terms of the power losses on the scale I = 103. We observe 
that the amplitude degrades faster with the density, i.e. for the same masked area, the loss 
of power is greater when more disks are placed than when the area of the disks is increased. 
The two cases have different geometries, in the first case the disks cover the map more 
uniformly, while in the second case larger contiguous unmasked areas are possible. 

Besides the change in amplitude in an individual mode, we also compare the two analyses 
in terms of their impact on the shape of the power spectrum. For this we consider two cases 
where the loss of amplitude in QL 103 amounts to about 50%. These are: a) n s = 0.60 deg -2 
for the fiducial 9 S , corresponding to a masked area / mas ked ~ 0.2 of the map, and b) 9 S = 
30 6> fwhm for the fiducial n s , corresponding to / mas kod ~ 0.3. The plot in Fig. [9] (left panel) 
shows that the loss of power is roughly scale-independent. In contrast, the variation of the 
size of the disks tilts the power spectrum shifting power to large scales (Fig. [9j right panel). 
The scales where the power increases are larger than the size of the enlarged reconstructed 
disk discussed in Sec. Ill, which for this example corresponds to i ~ 40. The shift of power 
is observed until 9 S ~ 40 #f w hm, beyond which the overall decrease in power overwhelms this 
tilt. 
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Figure 9: Shape degradation. Left panel: Ratio between the power spectra with and without 
masking reconstructed with the real-space (black line) and harmonic-space (gray line) estimators, 
for a mask of fiducial size and n s = 0.6 deg -2 density. Right panel: The same for a mask of fiducial 
density and size 0s = 30 6>f w h m - 



We also applied the harmonic-space estimator to these tests, even though it is outside 
the goal of our analysis to proper calibrate this estimator in the presence of masks. We 
found an increase of power of one to three orders of magnitude throughout all angular scales 
and spurious oscillations (Fig. [9]). When the masked area reaches 7% (26%) for the case of 
varying n s (6* s ), the power spectrum flattens and its amplitude starts to decrease with the 
masked area, throughout all scales. 



V. SUMMARY 

In this manuscript we study the effect of point source masks in the reconstruction of the 
weak lensing convergence from CMB maps. We used a variant of the minimum variance 
estimator with truncated support in real space. For comparison, we also used the original 
minimum variance estimator with full support in harmonic space. 

First we compared the reconstructions of convergence maps and power spectra made in 
complete versus masked temperature maps. The number of excisions in the masked map 
was determined by the number density of point sources reported in the literature, while the 
diameter of the masking disks was suggested by the size adopted in the ACT survey for 
modelling point sources. 

When using the real-space estimator, we found no difference in the quality of the esti- 
mation in the presence of masking. The input amplitude and shape were recovered after 
subtracting off the additive bias and applying a calibration vector derived from the recon- 
struction made in the absence of noise and masking. 

Conversely, when using the harmonic-space estimator, we found a decrease in the quality 
of the estimation of the convergence of around 40% at the map reconstruction level. We 
found that the amplitude and shape of the input power spectrum could not be recovered in 
the presence of noise or masking without further treatment of the calibration, since the simple 
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method of applying a calibration vector derived from the reconstruction in the absence of 
noise and masking cannot work in the presence of mode-coupling. There are several methods 
in the literature to deal with the problem of missing data. For example in Ref. [21] a method 
was developed to remove masking effects using a mode-coupling matrix. The "full" power 
spectrum is then calculated from the reconstructed "pseudo" power spectrum by means of 
a deconvolution involving matrix inversion, or with Bayesian methods. This was used in 
the recent first direct detection of a lensing potential in CMB data [3j. Another method 
consists in pre-whitening the map by applying an inverse covariance matrix weighting, in 
order to have a map that is less affected by mode-coupling. Such method was used in the 
first indirect detection of gravitational lensing in the CMB pp. 

We then proceeded to analyse the degradation of the reconstruction of the convergence 
power spectrum as a function of the parameters of the mask, varying the masked area 
between and 100%. 

We observe that the real space estimator yields less power as more map area is masked, 
with the power spectrum being typically recovered at 80% when the masked area reaches 
10% of the map. As a subdominant effect, we observe a shift of power from small to large 
scales, which is a consequence of the smearing of the masked disks observed in the histogram 
of the difference map. 

These observations strongly favour the use of the real-space estimator for the estimation 
of the weak lensing convergence in the case of maps with excision of points. The robustness 
of the real-space estimator, whether we vary the number density or the size of the masking 
disks, suggests that it can be used reliably on maps where extended areas have been masked. 

The computational cost of the real-space implementation scales with the number of pixels 
iV 2 on the CMB map and inversely with the pixel size #f w hm, whereas that of the harmonic- 
space implementation scales with N\og[N}. In particular, in the fiducial experiment used 
here, where N = 512 pixels and flfw h m = 7.2arcmin, the real-space estimator takes twice as 
long as our implementation of the harmonic-space estimator. Doubling the resolution, the 
real-space estimator takes four times longer. 

The problem of missing data is particularly relevant for the extraction of the polarization 
signal from the CMB, where it prevents a perfect E and B component separation [25J. 
Further studies can be envisaged to test the potentialities of this approach for lensing with 
CMB polarization. 
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